Evidence for universal four-body states tied to an Efimov trimer 
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We report on the measurement of four-body recombination rate coefficients in an atomic gas. Our results 
obtained with an ultracold sample of cesium atoms at negative scattering lengths show a resonant enhancement 
of losses and provide strong evidence for the existence of a pair of four-body states, which is strictly connected 
to Efimov trimers via universal relations. Our findings confirm recent theoretical predictions and demonstrate 
the enrichment of the Efimov scenario when a fourth particle is added to the generic three-body problem. 
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Few-body physics produces bizarre and counterintuitive 
phenomena, with the Efimov effect representing the major 
paradigm of the field |H]. Early in the 1970s, Efimov found 
a solution to the quantum three-body problem, predicting the 
existence of an infinite series of universal weakly bound three- 
body states. Surprisingly, these Efimov trimers can even exist 
under conditions where a weakly bound dimer state is absent 
J3. H 01- An essential prerequisite for the Efimov effect is 
a large two-body scattering length a, far exceeding the char- 
acteristic range of the interaction potential. Ultracold atomic 
systems with tunable interactions 15[] have opened up unprece- 
dented possibilities to explore such few-body quantum sys- 
tems under well controllable experimental conditions. In par- 
ticular, a can be made much larger than the van der Waals 
length r v dw the range of the interatomic interaction. 

In the last few years, signatures of Efimov states have been 
observed in ultracold atomic and molecular gases of cesium 
and recently in three-component Fermi gases of 
in a Bose gas of 39 K atoms 1 1 ill , and in mixtures 
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of 41 K and S/ Rb atoms 111211 . In all these experiments, Efimov 
states manifest themselves as resonantly enhanced losses, ei- 
ther in atomic three-body recombination or in atom-dimer re- 
laxation processes. The recent observations highlight the uni- 
versal character of Efimov states, and they also point to a rich 
playground for future experiments. 

As a next step in complexity, a system of four identical 
bosons with resonant two-body interaction challenges our un- 
derstanding of few-body physics. The extension of univer- 
sality to four-body systems has been attracting increasing in- 
terest both in theory H [3 H [H H Qj and experi- 
ment II 1 311 . A particular question under debate is the pos- 
sible relation between universal three- and four-body states 



FIG. 1 : (color online) Extended Efimov scenario describing a univer- 
sal system of four identical bosons; Energies are plotted as a function 
of the inverse scattering length. The red solid lines illustrate the pairs 
of universal tetramer states (Tetral and Tetra2) associated with each 
Efimov trimer (T). For illustrative purposes, we have artificially re- 
duced the universal Efimov scaling factor from 22.7 to about 2. The 
shaded regions indicate the scattering continuum associated with the 
relevant dissociation threshold. The four-body threshold is at zero 
energy and refers to four free atoms (A+A+A+A). In the a > re- 
gion, the dimer-atom-atom threshold (D+A+A) and the dimer-dimer 
threshold (D+D) are also depicted. The weakly bound dimer, only 
existing for a 2> r vc jw > 0, has universal halo character and its bind- 
ing energy is given by h 2 /(ma 2 ) B fTHl . The open arrow marks the 
intersection of the first Efimov trimer (T) with the atomic threshold, 
while the filled arrows indicate the corresponding locations of the 
two universal tetramer states. 
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In this context, Hammer and Platter 
predicted the four-body system to support universal tetramer 
states in close connection with Efimov trimers fl7Tl . 

Recently, von Stecher, D'Incao, and Greene presented key 
predictions for universal four-body states lfl9ll . For each Efi- 
mov trimer, they demonstrate the existence of a pair of uni- 
versal tetramer states according to the conjecture of Ref. lfl7ll . 
Such tetramer states are tied to the corresponding trimer 
through simple universal relations that do not invoke any four- 



body parameter 1 14, T(| IH]. The authors of Ref. lfl9ll suggest 
resonantly enhanced four-body recombination in an atomic 
gas as a probe for such universal tetramer states. They also 
find hints on the existence of one of the predicted four-body 
resonances by reinterpreting our earlier recombination mea- 
surements on 133 Cs atoms at large negative scattering lengths 
J3l. In this Letter, we present new measurements on the Cs 
system dedicated to four-body recombination in the particu- 
lar region of interest near a triatomic Efimov resonance. Our 
results clearly verify the central predictions of Ref. S. We 
observe two loss resonances as a signature of the predicted 
tetramer pair and we find strong evidence for the four-body 
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FIG. 2: (color online) Recombination losses in an ultracold sample 
of Cs atoms, (a) Loss fraction for a 50-nK sample after a storage 
time of 250 ms. Here we present all individual measurements to give 
an impression of the scatter of our data.The broad maximum at about 
— 870«o is caused by a triatomic Efimov resonance Q] and the shaded 
area highlights the resonant loss enhancement that we attribute to the 
four-body state Tetral. The three very narrow loss features (open cir- 
cles) are caused by known g-wave Feshbach resonances (20ll . which 
are irrelevant in the present context, (b) Loss fraction for a 30-nK 
sample after a storage time of 8 ms. Each data point represents 
the average values resulting from five individual measurements for 
a given a together with their statistical errors. The loss enhancement 
at around — 730an is caused by the state Tetra2. The solid lines are 
spline interpolations guiding the eye. 



nature of the underlying recombination process. 

The four-body extended Efimov scenario |T7, fl9Tl is 
schematically illustrated in Fig.Q] where the tetramer states 



(Tetral and Tetra2) and the relevant thresholds are depicted 
as a function of the inverse scattering length 1 /a. Within 
the four-body scenario, the Efimov trimers (T) are associated 
with trimer-atom thresholds (T+A, dashed lines). The pair 
of universal tetramer states (solid lines) lies below the corre- 
sponding T+A threshold. The four-body breakup threshold 
(A+A+A+A) defines zero energy and refers to the continuum 
of four free atoms. For completeness, we also show the a > 
region. Here, the picture is even richer because of the presence 
of the weakly bound dimer state, whichs leads to the dimer- 
atom-atom threshold (D+A+A) and the dimer-dimer threshold 
(D+D). In the four-body scenario, the tetramer states emerge 
at the atomic threshold for a < and connect to the D+D 
threshold for a > 0. 

The Efimov trimer intersects the atomic threshold at a = 
a j, which leads to the observed triatomic resonance J3]. The 
corresponding tetramer states are predicted lfl9ll to intersect 
the atomic threshold at scattering length values 



«Tetrai ~ 0.43 a\ and 4 etra2 « 0.9 a\. 



(1) 



These universal relations, linking three- and four-body reso- 
nances, express the fact that no additional parameter, namely 
the so-called four-body parameter, is needed to describe the 
system behavior. In contrast to the connection between uni- 



versal two- and three-body systems, where a three-body pa- 
rameter is required to locate the trimer states, the universal 
properties of the four-body system are thus directly related to 
the corresponding three-body subsystem. 

In analogy to the well-established fact that Efimov trimers 
lead to loss resonances in an atomic gas |0, Eltl . univer- 
sal four-body states can also be expected to manifest them- 
selves in a resonant increase of atomic losses 11911 . Resonant 
coupling between four colliding atoms and a tetramer state 
(a ~ fljetra < 0) drastically enhances four-body recombination 
to lower lying channels. Possible decay channels are trimer- 
atom, dimer-dimer, and dimer-atom-atom channels. In each 
of these recombination processes, we expect all the particles 
to rapidly escape from the trap, as the kinetic energy gained 
usually exceeds the trap depth. 

We prepare an ultracold optically trapped atomic sample 
in the lowest hyperfine sublevel (F = 3,mp = 3) lf22Tl . as de- 
scribed in Ref. 1 1311 . By varying the magnetic field between 6 
and 17 G, the scattering length a can be tuned from —1 100 to 
Oao 12311 . where ciq is Bohr's radius. For presenting our ex- 
perimental data in the following, we convert the applied mag- 
netic field into a using the fit formula of Ref. 0] • After several 
cooling and trapping stages 12311 . the atoms are loaded into an 
optical trap, formed by crossing two infrared laser beams 1 1311 . 
The trap frequencies in the three spatial directions are about 
{(O x ,(Oy,(O z ) =2tzx (10,46,65) Hz. Similar to H, we sup- 
port the optical trap by employing a magnetic levitation field 
acting against gravity. Evaporative cooling in the levitated 
trap is stopped just before the onset of Bose-Einstein conden- 
sation in order to avoid implosion of the gas. For our typical 
temperature of 50 nK, we obtain about 8 x 10 4 non-condensed 
atoms with a peak density of about 7 x 10 12 cm -3 . 

In a first set of experiments, we record the atom number 
after a fixed storage time in the optical trap for variable scat- 
tering length in the a < region. Figure|2]shows the observed 
losses, containing both three- and four-body contributions. 
The three-body part consists of a background that follows a 
general a 4 -scaling behavior lull l24l 12511 and resonant losses 
caused by the triatomic Efimov resonance, which for a 50-nK 
sample was observed to occur at a* T = — 870(10)ao 10 l26ll : 
this is consistent with the large losses shown in Fig.|2ja). Be- 
side this expected behavior of the three-body subsystem, we 
clearly observe two additional loss features, one located at 
about — 410flo [Fig.|2ja)] and one at about — 730«o [Fig. 12b)]. 
The observation of the resonance at — 730«o is particularly 
demanding and requires a careful choice of parameters as the 
signal needs to be discriminated against the very strong back- 
ground that is caused by three-body losses. Here we use a 
much shorter hold time of 8 ms, which is the shortest possible 
time required to ensure precise magnetic field control in our 
apparatus. 

We interpret the two observed resonant loss features as the 
predicted pair of four-body resonances |fl9ll . For the reso- 
nance positions we find antral / fl T — 0-47 and flxeti- a 2/ fl T — 
0.84, which are remarkably close to the predictions of Eq. (Q]). 

In a second set of experiments, we study the time- 
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FIG. 3: (color online) Time evolution of the number of atoms in 
an optically trapped sample. The solid lines are the fit to the data 
based on the full numerical solution of Eq. l[2}. The dashed and dot- 
ted lines correspond to a pure three-body and pure four-body de- 
cay, respectively (see text), (a) For dominant three-body collisions 
(a = — 870aoX 0>) an intermediate situation (a = — 510aoX ar >d (c) 
dominant four-body collisions (a = — 410ao)- 



dependence of the atomic decay in the optical trap. Here we 
focus on the region around the resonance at antral — — 410ao, 
where the three-body losses are comparatively weak and thus 
allow for a detailed analysis of the loss curves. Representative 
loss measurements for three different values of a are shown in 
Fig.S 

The observed decay can be fully attributed to three-body 
and four-body recombination collisions. This is due to the 
fact that inelastic two-body collisions of atoms in the lowest 
Zeeman sub-level are energetically suppressed, and one-body 
losses, such as background collisions or light-induced losses, 
can be completely neglected under our experimental condi- 
tions. The corresponding differential equation for the decay- 
ing atom number reads as 



N/N = -L 3 (n 2 )-L 4 (n 3 



(2) 



where L3 and L4 denote the three- and the four-body re- 
combination rate coefficient, respectively. The average den- 
sity is calculated by integrating the density over the volume 
(n 2 ) = {l/N)Jn 3 d 3 r and (« 3 ) = (l/N)Jn 4 d 3 r. By consid- 
ering a thermal density distribution of gaussian shape in the 



three-dimensional harmonic trap, we obtain (n ) = n 2 /\/T7 
and (« 3 ) = n 3 /8, with n p = \/%N[mW 2 /{Anks,T)fl 2 the peak 
density. Here, m is the atomic mass, T the temperature, and 
W = (OvfflyO)-) 1 / 3 the mean trap frequency. We determine the 
trap frequencies and the temperature by sloshing mode and 
time-of-fiight measurements, respectively. 

In general Eq. (0 is not analytically solvable. An analytic 
solution can be found in the limit of either pure three-body 
losses or pure four-body losses. 

Therefore we fit our decay curves with a numerical solution 
of Eq. (0, keeping both L3 and L4 as free parameters. Note 
that we have not included anti-evaporation heating 12511 in our 
model because we do not observe the corresponding tempera- 
ture increase in our experiments. We believe that, for the fast 
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FIG. 4: (color online) Loss rate coefficient for (a) three- and (b) four- 
body recombination as a function of the scattering length a. The 
measurements are taken at temperatures of about 40 nK. The values 
are obtained by fitting the numerical solution of Eq. (f2]l to the decay 
curve. The error bars on L3 and L4 are the statistical uncertainties 
from the fit evaluated with a resampling method (2711 . The open cir- 
cles in (a) refer to previous data at 250 nKp. The solid curves result 
from the theoretical model of Refs. IT^ . l28tl . 



decay observed here, the sample may not have enough time to 
thermalize. 

Our experimental data clearly reveal a qualitative change 
of the decay curves when a is tuned between a T and a Xetral . 
Fig. [3 a) shows that for a sa a T the loss is dominated by three- 
body recombination; here the full numerical fitting curve fol- 
lows the pure three-body solution. A different situation is 
found at — 410«o; see Fig-B c )- Here a pure three-body analy- 
sis cannot properly describe the observed behavior and the full 
numerical solution reveals a predominant four-body character. 
In intermediate situations, for which an example is shown in 
Fig. Ob), both three- and four-body processes significantly 
contribute to the observed decay. 

From the decay curves taken at different values of a we de- 
termine L3 and L4; the results are shown in Fig.HJa) and (b), 
respectively. The three-body contribution L3 follows previ- 
ously observed behavior J3l, as dictated by the a 4 -scaling in 
combination with the Efimov effect. 

Our major result is shown in Fig.@tb), where we plot the 
rate coefficient L4. Our data provide the first available quan- 
titative information on L4, establishing the role of four-body 
collisions in ultracold gases. For \a\ < |a Tetral l> wnere no uni- 
versal tetramer states exist, the four-body losses are typically 
very weak. Here, we measure L4 ~ 0.2 x 10~ 37 cm 9 /s. With 
increasing \a\, the system undergoes a significant change in 
its behavior, with four-body collisions dominating the atomic 
decay; see Fig. [2c). We observe a sharp increase of L4, which 
reaches its maximum value at a = — 412(2)«o- This observa- 
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tion is another strong piece of evidence for the predicted uni- 
versal four-body state at flj etral 01911 . To directly estimate the 
relative contributions of three- and four-body recombination, 
one can compare L3 with «o^4> where hq ~ 1.0 x 10 13 cm~ 3 
is the initial peak density at 40 nK. At resonance, «o^4 ex_ 
ceeds L3 by more than one order of magnitude, with L3 = 
0.7 x 10- 25 cm 6 /s and n L 4 = 3 x 10- 24 cm 6 /s. With fur- 
ther increasing \a\, L4 decreases and L3 increases such that 
L3 > /10L4. For \a\ > 700ao, the very fast three-body decay 
renders the analysis of the loss curves in terms of L4 unreli- 
able. In addition to the large statistical fit errors seen in our 
data in this region, other systematic error sources like a non- 
thermal evolution of the atomic density distribution can have 
a strong influence. 

Figured] also includes the theoretical predictions for L3 and 
L4 at 40 nK and demonstrates a remarkable qualitative agree- 
ment with our experimental results. The theoretical approach 
utilizes a solution of the four-body problem in the hyperspher- 
ical adiabatic representation lfl9ll : the derivation and associ- 
ated calculations of L4, adapted from lEitl . provide the first 
quantitative description of the four-body recombination rate. 
The calculations only require to fix the position of the tri- 
atomic Efimov resonance as determined in the previous ex- 
periment at 10 nK of Ref. |7J. The difference in the width and 
the amplitude of the four-body resonance between experimen- 
tal and theoretical data may be explained by different coupling 
to possible decay channels. 

Our work leads to important conclusions related to the con- 
cept of universality with increasing complexity. The observa- 
tion of the two four-body resonances close to the predicted po- 
sitions 1 19j] points to the universal character of the underlying 
states. This also supports the view or Refs. Q [H, O that 
a four-body parameter is not required to describe the system. 
Universal four-body states then emerge as a genuine conse- 
quence of the Efimov spectrum. This also provides a novel 
way to test Efimov physics. The Efimovian character of a 
three-body resonance can be probed by observing the univer- 
sal tetramer resonances tied to it, without the necessity to ex- 
plore the full geometric scaling of Efimov physics by chang- 
ing the scattering length by orders of magnitude. 

While our present work has focussed on four-body phe- 
nomena at negative scattering length, a further exciting step 
will be the exploration of the entire four-body spectrum. For 
positive scattering lengths, the spectrum becomes richer and 
new phenomena can be expected such as resonant interac- 
tions between four-body states and two-dimer states. In this 
way, experiments on few-body phenomena in ultracold atoms 
will keep on challenging our understanding of the universal 
physics of a few resonantly interacting particles. 

We are aware of related results in 39 K, in which enhanced 
losses near a Feshbach resonance may be interpreted as a four- 
body resonance l29tl . 
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